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Abstract 
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1 Introduction 



The phase diagrams of many physical systems have a "critical region" where all traditional 
approximation methods, such as mean-field theory and its generalizations, fail completely 
to provide an accurate description of the system's behavior. This is due to the existence of 
a "critical point," approaching which, some statistical-mechanical quantities diverge, while 
others stay finite but have divergent derivatives. Experimentally, it is found that these 
quantities usually behave in the critical region as a power law. The exponents that appear 
in those power laws, called critical exponents, describe the nature of the singularities 
at the critical point. The theory of critical phenomena based on the renormalization 
group suggests that statistical-mechanical systems fall into "universality classes" such 
that systems belonging to the same universality class have the same critical exponents. 

It should be emphasized that the phenomenon of universality was discovered at least 
twice prior to the introduction of the renormalization group. In fact, in what Alan 
Sokal [IS] calls the Dark Ages of the theory of critical phenomena, before the 1940 's, 
physicists generally believed that all systems have the same critical exponents, namely 
those of the mean-field theory of Weiss [17] (1907) or its analogue for fluids, the van der 
Waals theory (1873). The failure of this type of universality became apparent as 
early as 1900, following experiments on fluid systems, but it began to be taken seriously 
only after Onsager's exact solution [37] (explicitly displaying non-mean-field exponents) 
in 1944 and the rediscovery of the experimental evidence of non-mean-field values for 
critical exponents by Guggenheim in 1945 JH]- Nonetheless, although mean-field theory 
is clearly incorrect for short range models, some universality does seem to hold. Many, if 
not all, different fluids, for instance, seem to have the same value for the critical exponent 
(3 (related to density fluctuations), and it is believed that, for example, carbon dioxide, 
xenon and the three dimensional Ising model should all have the same critical exponents. 
Maybe even more surprisingly, it was soon realized that some (but not all) magnetic sys- 
tems have the same critical exponents as do the fluids. This remarkable phenomenon 
seems to suggest the existence of a mechanism that makes the details of the interaction 
irrelevant in the critical region. Nevertheless, the critical exponents should depend on 
the dimensionality of the system and on any symmetries in the Hamiltonian. Despite 
being a very plausible and appealing heuristic idea, backed up by renormalization-group 
arguments and empirical evidence, only very few cases are known in which universality 
has actually been proved. 

Percolation, with its simplicity and important physical applications, is a natural can- 
didate for studying universality. This is especially so after the ground-breaking work of 
Lawler, Schramm, and Werner [23IHI23I13I27JI2HI1EIII1, who identified the only 
possible conformally invariant scaling limit of critical percolation and derived many of its 
properties, and that of Smirnov [121113, who proved that, indeed, critical site percolation 
on the triangular lattice has a conformally invariant scaling limit. The combination of 
those results made it possible [Hj to verify the values of the critical exponents predicted 
in the physical literature in the case of critical site percolation on the triangular lattice 
(and to derive also some results that had not appeared in the physics literature, such 
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as an analogue of Cardy's formula "in the bulk" |41j, or the description of the so-called 
backbone exponent [27]). 

It is generally accepted that the lattice should play no role in the scaling limit, and 
that there should be no difference between bond and site percolation. In other words, 
two-dimensional critical (independent) percolation models, both site and bond, should 
belong to the same universality class, regardless of the lattice (at least for periodic lattices 
like the square, triangular or hexagonal lattice). Once again, though, despite being a 
very natural and plausible conjecture, such universality has not yet been proved. There 
is however another natural direction in which to study universality, which consists in 
analyzing critical percolation models on a given lattice that differ in their dependence 
structures. It is this direction that we pursue in this paper. 

The cellular automata that we use to generate our families of dependent percolation 
models arise naturally in the study of the zero-temperature limit of Glauber dynamics or 
as coarsening or agreement-inducing dynamics. The action of such cellular automata can 
be viewed as a sort of "small (local) perturbation" of the original, independent percolation 
model, and our main corollary can be viewed as proving a form of universality for two 
dimensional percolation. Therefore, we provide an explicit example of the principle of 
universality, in the strong form concerning scaling limits. 

To be more precise, there are at least two, a priori different, notions of universality, 
one concerning the critical exponents discussed above, and a second one concerning the 
continuum scaling limit. The two concepts are closely related, but in this paper we are 
only concerned with the second type of universality (the first type will be considered in a 
future paper). 

We consider a family of dependent percolation models that arise through a (discrete 
time) deterministic cellular automaton T acting on site percolation configurations a on the 
set of sites of the triangular lattice T. Each configuration a corresponds to an assignment 
of —1 or +1 to the vertices of T. The variable o~ x , corresponding to the value of a at 
x, is commonly called a spin variable. At discrete times n = 1,2, . . ., each spin a x is 
updated according to the following rules (later, in Section 01 we will introduce other 
cellular automata, both on T and H, generated by different rules): 

• if x has three or more neighbors whose spin is the same as a x , then the latter does 
not change value, 

• if x has only two neighbors y\ and ?/2 that agree with x, and y\ and y<i are not 
neighbors, then o x does not change value, 

• otherwise, a x changes value: a x — > —cr x . 

The starting configuration a , at time n — 0, of our cellular automaton is chosen from 
a Bernoulli product measure corresponding to independent critical percolation, and the 
distributions at times n > 1 (including the final state as n — ► oo) of the discrete time 
deterministic dynamical process a n are the other members of our family of dependent 
percolation models. We show that all those dependent percolation models have the same 
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scaling limit, thus providing an explicit example in which universality can be proven. This 
comes as a corollary of our main result, Theorem ^ To explain the main result, we first 
need some terminology 

In the scaling limit, the microscopic scale of the system (i.e., the lattice spacing 5) is 
sent to zero, while focus is kept on features manifested on a macroscopic scale. In the 
case of percolation, it is far from obvious how to describe such a limit and we only do so 
briefly here; for more details, see [TJI21I3]. We will make use of the approach introduced by 
Aizenman and Burchard [3] (see also jl]) applied to portions of the boundaries between 
clusters of opposite sign, and present our results in terms of closed collections of curves 
in the one-point compactification M 2 of M 2 , which we identify (via the stereographic 
projection) with the two-dimensional unit sphere. For each fixed 5 > 0, the curves are, 
before compactification, polygonal paths of step size 5 (i.e., polygonal paths between 
sites of dual lattice). The distance between curves is defined so that two curves are 
close if they shadow each other in a metric which shrinks at infinity (for the details, see 
Sections O and 

The dynamics allows one to construct the whole family of percolation models for all 
n G {0, 1, ... , oo} on the same probability space, i.e., there is a natural coupling u, realized 
through the dynamics, between any two percolation models in the family. In terms of this 
coupling, Theorem Q states, roughly speaking, that the //-probability that the distance 
between the two collections of curves corresponding to two distinct percolation models in 
the family is bounded away from zero vanishes as 5 — > 0. This means that, in the limit 
5 — > 0, given any two models in the family, for every curve in one of them, there exists a 
curve in the other one that shadows the first curve and vice-versa. 

The rest of the paper is organized as follows. In Section we describe the behavior of 
the cellular automaton T, define the family of dependent percolation models dynamically 
generated by T, and state the main results. The proofs of the main results are contained 
in Section El The dynamics described in Section |21 is chosen as a prototypical example, 
but is not the only one for which our results apply. In Section 01 we introduce other such 
dynamics (both on T and H), which can be obtained as suitable zero-temperature limits 
of stochastic Ising models. 

2 Definitions and results 

We start this section by giving a more detailed definition of one of the families of dependent 
percolation models that are the object of investigation of this paper. The models in that 
family are defined on the triangular lattice T, embedded in M 2 by identifying the sites 
of T with the elementary cells (i.e., regular hexagons) of the hexagonal (or honeycomb) 
lattice H (see Figured)). We will use those models as a paradigm and will give for them 
explicit and detailed proofs of the results. Later on, we will point out how to modify the 
proofs to adapt them to the other models discussed in the paper. 

In the rest of the paper, points of M 2 will be denoted by u and v, while for the sites 
of T we will use the Latin letters x, y, z and the Greek letters ( and £. An edge of T 
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incident on sites x and y will be denoted by r} X:V , while by if x y we denote the dual edge 
(in H) perpendicular bisector of r} x>y . 




Figure 1: The hexagonal (or honeycomb) lattice. 



2.1 The dynamics and the percolation models 

We construct a family of dependent percolation models by means of a cellular automaton 
T acting on site percolation configurations on the triangular lattice T, i.e., T : Q — > Q, 
where Q is the set of configurations {— 1,+1} T . This family is parametrized by n G 
{0, 1, . . .} U {oo} representing (discrete) time. The initial configuration cr° consists of an 
assignment of — 1 or +1 to the sites of T. At times, we will identify the spin variable a x with 
the corresponding site (or hexagon) x. We choose a° according to a probability measure /j, 
corresponding to independent identically distributed <r°'s with /z°(cr)j! = +1) = A G [0, 1]. 
With the exception of Proposition 12.21 we will set A = 1/2, so that /j, is the distribution 
corresponding to critical independent site percolation. We denote by \i n the distribution 
of a n and write fi n = T fi°, where T is a map in the space of measures on Q. The action 
of the cellular automaton T can be described as follows: 

{a™ if x has at least two neighbors y± and yi such that cr™ = a™ 2 = a™, 
and y\ and y2 are not neighbors 
— a™ otherwise 

(1) 

Once the initial percolation configuration <r° is chosen, the dynamics is completely deter- 
ministic, that is, T is a deterministic cellular automaton with random initial state (see 
Figs. El and EJ). 

Certain configurations are stable for the dynamics, in other words, they are absorbing 
states for the cellular automaton. To see this, let us consider a loop in T expressed 
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as a sequence of sites (Co, • • • , Cfc) which are distinct except that Co = Cfc an d suppose 
moreover that k > 6 and that and Q + i are not neighbors; we will call such a sequence 
an m-loop. If a^ — — . . . — cr^ k _ 1 then every site Q in such an m-loop has two 
neighbors, and Q+i, such that and are not neighbors of each other and 
cr^_ a = cr^ = cr^ i+1 . According to the rules of the dynamics, is therefore stable, that 
is, retains the same sign at all future times. Other stable configurations are "barbells," 
where a barbell consists of two disjoint such m-loops connected by a stable m-path. The 
stability of certain loops under the action of T will be a key ingredient in the proof of the 
main theorem. A more precise definition of m-paths and m-loops is given in Section 01 



Figure 2: Example of local configuration with unstable spins. The numbered hexagons 
correspond to spins that will flip and the numbers indicate at what time step the spin 
flips occur. 



An important feature of this dynamics is that almost surely every spin flips only a 
finite number of times and every local configuration gets fixated in finite time. To show 
this, we introduce a formal Hamiltonian 



(x,y) x 

where Y^,( x y ) denotes the sum over all pairs of neighbor sites, each pair counted once, and 



T~tx{v) = - CfxOy, (3) 

yeN(x) 

where J\f(x) is the set of six (nearest) neighbors of x. We also introduce a "local energy" 

(x,y) zedA x£A 

x,yeA x£Af(z) 



Figure 3: Same local configuration after all the unstable spins have flipped. 



where A is a subset of T and <9A is the outer boundary of A, i.e., {( ^ A : x G 
J\f(C) for some x G A}. Notice that although the total energy TC(a) is almost surely 
infinite and is therefore only defined formally, we will only use local energies of finite 
subsets of T. 

The notion of the energy change caused by a spin flip is somewhat ambiguous in a 
cellular automaton because of the synchronous dynamics and hence multiple simultaneous 
spin flips. Nevertheless, it is easy to show (see the proof of Proposition [OJ that each step 
of the dynamics either lowers or leaves unchanged the energy - both locally and globally. 
In this sense, our cellular automaton can be considered a zero-temperature dynamics (see, 
for example, [32] and references therein). 

Proposition 2.1. For all values of X, almost surely, every spin flips only a finite number 
of times. 

Proof. By the translation invariance and ergodicity of the model, it is enough to prove 
the claim for the origin. At time zero, and for all values of A, the origin is almost surely 
surrounded by an m-loop T of spins of constant sign. (For A = 1/2, there are infinitely 
many loops of both signs surrounding the origin.) Such an m-loop is stable for the 
dynamics and its spins retain the same sign at all times. Call A the (a.s. finite) region 
surrounded by T. The energy H\ of the finite region A is bounded below. Consequently, 
if we can show that no step of the dynamics ever raises 7i\, it would follow that there 
can only be a finite number of steps that strictly lower the energy 7i\. 

Call an edge r} x , y satisfied if a x = a y and unsatisfied otherwise; then the change in 
energy TC\{o' n+1 ) — 7Ya(o" ?1 ) is twice the difference between the number of satisfied edges at 
time n that become unsatisfied at time n+1 and the number of unsatisfied edges at time 
n that become satisfied at time n + The edges that change from satisfied to unsatisfied 
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or vice- versa are those between spins that flip and their neighbors that do not, so 

waC^ 1 ) - w A (o = E E - = E A «^> ( 5 ) 

xeA : <t£+ V< 2/eAf(z) xeA : c% +1 

where 

A n H x = n x (<T n+1 )-n x ((T n ). (6) 

Notice that the only nonzero contributions in the first sum of (0) come from those sites 
y G A that do not flip at time n. We want to show that A n TC x < for all x G A and find 
some y with A n H y < (assuming there was at least one spin flip inside A at time n). 

Call D x the number of disagreeing neighbors of x at time n and notice that a necessary 
condition for the spin at site x to flip at that time is D™ > 4. Let us first consider the case 
D x > 5 and assume, without loss of generality, that a™ = — 1 and cr™ +1 = +1. Then, at 
time n, site x has at least five plus-neighbors, and at least three of them have plus-spins 
at time n + 1 (those having at time n two plus-neighbors that are not neighbors of each 
other). This implies that the number of edges incident on x that change from unsatisfied 
to satisfied is at least three, while the number of edges that change from satisfied to 
unsatisfied is one. Then, 

A n H x = E W ~ < 0. (7) 

Next, we consider the case D" = 4. Again, we can assume that a x = — 1 = — er" +1 . In 
this case, one can have two types of spin flips, one with TC x (a n+1 ) — H. x (cr n ) < and one 
with Hx(& n+1 ) — 7~L x (c> n ) = 0. The second type occurs when, of the four neighbors of x 
that are plus at time n, two remain plus at time n + 1 (we call them y\ and y-i) and the 
other two flip to minus (we call them z\ and Z2), while the two neighbors that are minus 
at time n remain minus at time n + 1 (we call them yi and y^). In this situation, site x 
disagrees with four of its six neighbors both at time n and at time n + 1 (see Figure 0J) 
and therefore H x (a n+1 ) - H x (a n ) = 0. 

The spins at Z\ and z 2 flip together with the spin at x. Each energy, Tt Zl and TC Z2 , is 
lowered or left unchanged. If either is lowered, then the energy 7i\ is lowered. If neither 
is lowered, then Z\ and z 2 must be in the same situation as x (but with minus and plus 
interchanged), which requires a configuration that looks locally like the one in the left 
part of Figure |U (or one equivalent to it under some lattice symmetry). However, such 
a local configuration cannot extend forever; it must be finite and contained in A. This 
implies that we will necessarily find at least one spin y (in fact, at least two spins) that 
flip together with the spin at x and such that A n H y < 0. 
Thus, if at time n some site x G A flips, 

H A (a n+1 )-H A ( ( T n )= E ^nH x <0. (8) 

It follows that there can only be a finite number of times n at which spins in A (in par- 
ticular, the origin) flip. This completes the proof. □ 
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Time n 



Time n+ 1 



Figure 4: Example of a step of the dynamics acting on a local configuration that leaves 
the energy of some sites (including x) unchanged and decreases that of other sites (the 
three spins at the top and the three at the bottom are stable, as determined also by spins 
that do not appear in the figure). 



Let us now give two results that are analogous to results proved in [5] for a related 
cellular automaton that will be discussed below in Section 0] 

Proposition 2.2. If X > 1/2 (respectively, < 1/2), then for almost every a , there is 
percolation of + 1 (respectively, —1 ) spins in a n for any n G [0, oo]. 

Proof. We only give the proof for A > 1/2, since the case A < 1/2 is the same by sym- 
metry. If A > 1/2, since the critical value for independent (Bernoulli) percolation on the 
triangular lattice is exactly 1/2, there is at time zero percolation of +1 spins. This implies 
the existence, at time 0, of doubly-infinite plus-paths that are stable. Therefore, there is 
percolation of +1 spins for all n > 0. □ 

We denote by C x the cluster at x (for a configuration a), i.e., the maximal connected 
set B G T such that x G B and a y = a x for all y G B. We write C x {n) to indicate the 
cluster at x for a n . C x {n) is random and its distribution is denoted fi n . E^n denotes 
expectation with respect to //". 

Proposition 2.3. For A = 1/2, the following two properties are valid. 

1. For almost every a , there is no percolation in a n of either +1 or — 1 spins for any 
n G [0, oo]. 

2. The mean cluster size in a n is infinite for any n G [0, oo] : for any x G T, 

E^(|C x .|)=E M o(|C x .WI) = oo. (9) 
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Proof. To prove the first claim, notice that at time zero the origin is almost surely 
surrounded by both a plus and a minus m-loop. Those loops are stable and prevent the 
cluster at the origin, be it a plus or a minus cluster, from percolating at all subsequent 
times. Therefore the probability that the origin belongs to an infinite cluster is zero for 
all n > 0. 

To prove the second claim, we note that because of the absence of percolation for 
either sign, it follows from a theorem of Russo [HH] applied to the triangular lattice, that 
the mean cluster size of both plus and minus clusters diverges. □ 

Before stating our main theorem, we need some more definitions to formulate the 
continuum scaling limit. We adopt the approach of |3] (see also jl]). 

2.2 Compact ificat ion of IR 2 

The scaling limit can be taken by focusing on fixed finite regions, A C K 2 , or by 

treating the whole IR 2 . The second option is more convenient, because it avoids technical 
issues that arise near the boundary of A. 

A convenient way of dealing with the whole IR 2 is to replace the Euclidean metric with 
a distance function d(-, •) defined on IR 2 x IR 2 by 

d{u, v) = inf J (1 + I0I 2 )- 1 ds, (10) 

where the infimum is over all smooth curves 0(s) joining u with v, parametrized by arc- 
length s, and | • | denotes the Euclidean norm. This metric is equivalent to the Euclidean 
metric in bounded regions, but it has the advantage of making IR 2 precompact. Adding a 
single point at infinity yields the compact space M 2 which is isometric, via stereographic 
projection, to the two-dimensional sphere. 

2.3 The space of curves 

Denote by S the complete, separable metric space of continuous curves in IR 2 with a 
distance D(-, •) based on the metric defined by eq. (fTUj) as follows. Curves are regarded 
as equivalence classes of continuous functions ^(t) from the unit interval to IR 2 , modulo 
monotonic reparametrizations. C will represent a particular curve and 7 (i),£ G [0,1], a 
particular parametrization of C, while JF will represent a set of curves. The distance D 
between two curves, C\ and C2, is defined by 

D(C a ,C 2 ) = inf sup d( 7l (/i(t)),72(/ 2 (t))), (11) 
a>/2 te[o,i] 

where 71 and 72 are particular parametrizations of C\ and C2, and the infimum is over the 
set of all monotone (increasing or decreasing) continuous functions from the unit interval 
onto itself. The distance between two closed sets of curves is defined by the induced 
Hausdorff metric as follows: 

dist(J r , T') < e <=> VC e 3 C G T' with D(C, C) < e, and vice-versa. (12) 
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For each fixed S > 0, the random curves that we consider are polygonal paths in 
the hexagonal lattice 8M., dual of 8T, consisting of connected portions of the boundaries 
between plus and minus clusters in ST. A subscript 5 may then be added to indicate that 
the curves correspond to a model with a "short distance cutoff." The probability measure 
Hg denotes the distribution of the random set of curves consisting of the polygonal 
paths on SW generated by T n acting on cr°. 

2.4 Main results 

Since the cellular automaton is deterministic, all the percolation models a n for n > 
are automatically coupled, once they are all constructed on the single probability space 
(Q, S,//°) on which <r° is defined. The following theorem is valid for any A, but A = 1/2 
is the only interesting case, therefore we restrict attention to it. 

Theorem 1. For X = 1/2, the Hausdorff distance between the system of random curves 
at time and the corresponding system of curves at time n goes to zero almost 
surely as 5 — > 0, for each n G [1, oo]; i.e., for n° -almost every a , 

lim dist(T 5 °, 7%) = 0, for any n G [1, oo]. (13) 

5^0 

The main application of the theorem is that the scaling limits of our family JF™ of perco- 
lation models, if they exist, must be the same for all n G [0, oo]: 

Corollary 2.1. Suppose that critical site percolation on the triangular lattice has a unique 
scaling limit in the Aizenman-Burchard sense J$j, i.e., !Fg converges in distribution as 
5 — > (for A = 1/2). Then, for every n G [l,oo], JF™ converges in distribution to the 
same limit. 

3 Proofs 

In this section, we give the proofs of the main results. We start by reminding the reader of 
the definitions of m-path and m-loop and by giving some new definitions and two lemmas 
which will be used later. 

Definition 3.1. A path T between x andy in ST (embedded in M 2 ) is an ordered sequence 
of sites (Co = x, . . . , (k = y) with Q ^ Q for i j and Q + i G A/"(Ci) ■ A loop is a sequence 
(Co, • • • , Cfc+i = Co) with k > 1 such that (Co, • ■ • , Cfe) and (Ci, . . . , Cfc+i) are paths. 

We call a path (Co, • • • , Cfe) (resp., a loop (Co, • • • , Cfc+i — Co) ) an m-path (resp., an 
m-loop^) ifd-i and Ci+i are not neighbors for i — 1, . . . , k — 1 (resp., for i — 1, . . . , k+ 1, 
where ( k+2 = (i). 

Notice that for every path T between x and y, there always exists at least one m-path F 
between the same sites. 
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Definition 3.2. A boundary path (la-pathj V* is an ordered sequence = T]t ■ ■ ■ ,Vk ~ 

T)l k £ h ) °f distinct dual edges such that either Q + i = Q and G A/"(£i) or = and 
Ci+i G Af(Ci), an d = -cr 6 f or all i — 0, . . . ,k. 

We call a maximal boundary path simply a boundary; it can either be a doubly-infinite 
path or a finite loop. 

b-paths represent the (random) curves introduced in the previous section for a fixed 
value of the short distance cutoff 5. A collection of such curves in the compact space M? 
is indicated by 

b-paths T* are parametrized by functions j(t), with t G [0, 1]. When we write that, 
for t in some interval [ti,t 2 ] (the interval could as well be open or half-open), j(t) G rf x y 
we mean that the parametrization ^(t) for t between t\ and t 2 is irrelevant and can be 
chosen, for example, so that |^^| is constant for t G [ii, £2]- The notation T*(u, v), where 
u and v can be dual sites or generic points of IR 2 fl T*, stands for the portion of T* between 
u and v . 

Definition 3.3. A stable loop / (for some a) is an m-loop ((0 = x, . . . ,(k = x) such 
that o~£ = er^ = . . . = a^ k . 

Definition 3.4. We say that a dual edge rj* is stable if x and y belong to stable loops 
of opposite sign. 

Lemma 3.1. In general, a constant sign m-path (( , . . . ,(k) is 'fixated" (i.e., retains that 
same sign in a n for all < n < 00) if Co and are fixated 

Proof. The claim of the lemma is straightforward. 

Lemma 3.2. For A = 1/2, there is a one to one mapping from boundaries in to 
"parent" boundaries in 

Proof. Let T* G be a boundary at time n and int(r*) be the set of sites of T (i.e., 
hexagons) "surrounded" by T*. Let C the (unique) constant sign cluster contained in 
int(r*) that has sites next to T*. Suppose that sites x and y in C do not change sign at 
time n + 1, then, at that time, they must belong to the same cluster C. (This means 
that a cluster cannot split in two or more pieces under the effect of the dynamics.) The 
reason is that x and y are connected by an m-path T at time n because they are in the 
same cluster, and since they have not flipped between time n and time n + 1, Lemma f3. II 
(or more accurately, the single time step analogue of Lemma fo.lj) implies that the sites in 
r have not flipped either, so that x and y are still connected at time n + 1 and therefore 
belong to the same cluster. 

If at least one site x in C does not change sign at time n + 1, we say that C has 
survived (and evolved into a new cluster C that contains x at time n + 1). From what 
we just said, the sites of C that retain the same sign form a unique cluster C We call 
C the parent cluster of C and the external boundary T* G of C the parent of the 
external boundary G of C (note that C is a.s. finite - see Proposition 12. 3|) . 
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To prove that there is a one to one mapping between boundaries in J 7 ^ + and those in J-] 1 , 
it remains to show that the parent boundary of each element of is unique. (This 

means that clusters cannot merge under the effect of the dynamics.) 

If this were not the case, then a cluster C could have two or more distinct parent 
clusters, C%, C 2 , . . . . Notice that each of the parent clusters at time n is surrounded by a 
constant sign m-loop Tj, i = 1,2,..., which is stable. Suppose that x x G C\ and x 2 G C 2 
retain their sign at time n + 1. Without loss of generality, we shall assume that this sign 
is plus. Assuming that C\ and C 2 are both parents of C, X\ and x 2 should both belong 
to C at time n + 1 and therefore be connected by a plus-path. But this contradicts the 
fact that, at time n, x\ G int(ri) and x 2 G int(r 2 ), that is, they belong to the interiors of 
two disjoint stable minus-loops. □ 

Notice that boundaries cannot be "created," but can "disappear," as complete clusters 
are "eaten" by the dynamics. 

3.1 Proof of Theorem [1] 

Let us start, for simplicity, with the case of a single b-path. For any e > 0, given a 
b-path Tq at time 0, parametrized by some j(t), we will find a path T* at time n with 
parametrization j'(t) such that, for S small enough, 

sup d( 7 (t), 7 '(t)) < e + 2S (14) 
te[o,i] 

and vice- versa (i.e., given T* and 7', we need can find Tq and 7 so that (JHJ), in which 
the dependence on the scale factor S has been suppressed, is valid). Later we will require 
that this holds simultaneously for all the curves in JF5 and J-'g, as required by eq. (fT3)). 

Let B 1 (R) be the ball of radius R, B 1 (R) = {u G R 2 : \u\ < R} in the Euclidean 
metric, and B 2 (R) = {u G M 2 : d(u) < R} the ball of radius R in the metric fllOl) . For 
a given e > 0, we divide R 2 into two regions: B 1 (6/e) and M 2 \ B 1 (6/e). We start by 
showing that, thanks to the choice of the metric (|10jh one only has to worry about curves 
(or polygonal paths) that intersect B 1 (6/e). In fact, the distance between any two points 
u,!)£l 2 \ 5 1 (6/e) satisfies the following bound 

POO 

d(u, v) < d(u, 00) + d(u, 00) < 2 / [1 + (s + 6/e) 2 ]" 1 rfs < e/3. (15) 

Thus, given any curve in contained completely in M 2 \i? 1 (6/£:), it can be approximated 
by any curve in also contained in R 2 \ B l (6/e), and viceversa. The existence of such 
curves in is clearly not a problem, since the region M 2 \ 5 1 (6/£) contains an infinite 
subset of ST and therefore there is zero probability that it doesn't contain any b-path at 
time zero. There is also zero probability that it contains no stable b-path at time zero, 
but any such b-path also belongs to T§. 

Before we can proceed, we need the following lemma, which is a consequence of the 
the fact that at time zero we are dealing with a Bernoulli product measure. In this lemma 
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(and elsewhere), the diameter diam(-) of a subset of M 2 is defined by using the Euclidean 
metric. 

Lemma 3.3. Let rj^ be any (deterministic) dual edge; then for some constant c > 0, 

/j° (3T* 3 t]q : diam(r*) > M and T* does not contain at least one stable edge) < e~ cM . 

(16) 

Proof. To prove the lemma, we partition the hexagonal lattice into regions Qi as in Fig- 
ure El We then do an algorithmic construction of T*, starting from percolation 
exploration process, but with the additional rule that, when the exploration process hits 
the boundary of any Qi for the first time, all the hexagons in Qi are checked next (accord- 
ing to some deterministic order). From every entrance point of Qi, there is a choice of 
the values of the spins of the outermost layer of Qi that forces T* to enter Qi. Therefore, 
when T* hits Qi, it always has a positive probability of entering the region. 




Figure 5: Elementary cell for the partition of EI used in Lemma [3. 31 Notice that the cell 
is made out of seven smaller cells, each of them formed by seven hexagons. 

We call Fi the event that a (dual) stable edge is found inside Qi, belonging to T*. It 
is easy to see that such an event has positive probability, bounded away from zero by a 
constant that does not depend on how the exploration process enters the region Qi. In 
fact, from any entrance point, there is clearly a choice of the values of the spins in Qi that 
forces T* to cut Qi in two symmetric parts, containing spins of opposite sign. 

Now, if diam(r*) > M, then T* must clearly visit at least O(M) different regions Qi. 
The conclusion of the proof should now be clear (see, for example, |15j). □ 
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With this lemma, we can now proceed to the proof of the theorem. As explained 
before, we restrict attention to paths that intersect B 1 (Q/e). Given a b-path Tq = (rj* = 
n Co fc • • • ' = V( k £ k ) m with parametrization 7(t), call uq the point in M 2 where the 
first (dual) edge t/q begins. The following algorithmic construction produces a sequence 
Uo, . . . , un of points in Tq. 

1. Start with u . 

2. Once uq, . . . , U{ have been constructed, if Ui G B 1 (6/e) take the ball B 1 (e/3) cen- 
tered at Ui and of radius e/3 and let Ui+i be the first intersection of Tq \ Tq(uq, u,/) 
with dB*. (e/3), if Ui ^ B 1 (6/e) take the ball B 2 .(e/3) centered at and of radius 
e/3 and let u i+ i be the first intersection of Tq \ Tq(uq, Uj) with dB 2 i .(e/3). 

3. Terminate when there is no next Mj. 

During the construction of the sequence u , . . . ,un, Tq is split in N + 1 pieces, the 
first N having diameter at least e/3. The construction also produces a sequence of balls 
B^ o (e/3), . . . , Bl^(e/3), with ji = 1 or 2. Notice that no two successive Mj's can lie outside 
of B 1 (6/e). In fact, if for some i, U; L lies outside of B 1 (6/e), u i+1 belongs to dB^.(e/3), 
which is contained inside B 1 (Q/e), due to the choice of the metric. Each Ui contained in 
B 1 (Q/e) lies on an edge of SM, but no more than one Ui can lie on the same edge since Tq 
is self- avoiding and cannot use the same edge or site more than once. Also, the number 
of Mj's lying outside of B 1 (6/e) cannot be larger than (one more than) the number of the 
UiS lying inside B 1 (6/e). Therefore, iV < const x (e5)~ 2 . 

For any two successive balls, B£(e/3) and B J u i +\(e/3), let Oi = B^(e/3 + 5) U 
Bu^\(e/3 + 5). Now assume that there exists a sequence fjQ : ...fj%_ 1 of stable (dual) 
edges of Tq, with 77* contained in T^Ui, u i+ i). T^fj* , f}* +1 ) is contained in Oi (for fixed 
e and small enough 5). Also contained in Oi are the two paths (i.e., paths in 5T which 
may be thought as sequences of hexagons) whose hexagons are next to T^fj*, fj* +1 ), one 
on each side. From those two paths one can extract two subsets that are m-paths which 
are stable since the first and last hexagons of each one of them is stable (such hexagons 
must be stable since they are next to fj* and fj* +1 ). The two m-paths so constructed 
constitute a "barrier" that limits the movements of the boundary, so that r* (fj*, fj* +1 ) is 
in fact confined to lie within those two m-paths and thus within Oj. To parametrize T*, 
we use any parametrization ^'(t) such that ^'(t) = ^(t) whenever j(t) e fj*. Using this 
parametrization and the previous fact, it is clear that the distance between Tl(fj*,fj* +1 ) 
and r* (fj*, fj* +1 ) does not exceed e + 25, the diameter of 0«. Therefore, conditioning on the 
existence of the above sequence fj*, . . . ,fj* N _ x of stable (dual) edges of Tq, we can conclude 
that 

sup d( 7 (*),Y(f))<e + 25. (17) 
te[o,i] 

It remains to prove the existence of the sequence t]q, ■ ■ ■ °f stable (dual) edges. 

To do that, let us call the event that Tq(ui, ■Uj+i) does not contain at least one stable 
edge, and let A = U^q 1 ^ be the event that at least one of the first N pieces of Tq does 
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not have any stable edge. Then, considering that the total number of edges contained in 
B 1 (6/e) is bounded by const x (e5)~ 2 and using Lemma EHfl we have 

(i° s (A) < (s5)~ 2 e~ c ' {e/6) (18) 

for some d > 0. Equation (fTSj) means that the probability of not finding at least one 
stable edge in each of the first iV pieces of Tq is very small and goes to 0, for fixed e, as 
5^0. This is enough to conclude that, with high probability (going to 1 as 5 — > 0), 
equation (fTTj) holds. 

This proves one direction of the claim, in the case of a single curve. To obtain the 
other direction, notice that a large b-path at time n is part of a complete boundary T* 
which must come from a line of "ancestors" (see Lemma f3.2|) that starts with some Tq at 
time 0. Therefore, one can apply the above arguments to Tq, provided that the latter is 
large enough. Although the proof of this last fact is very simple, it is convenient to state 
it as a separate lemma. 

Lemma 3.4. Set 5 = 1 for simplicity; then for any boundary T* at time n, there is an 
ancestor Tq, with diam{TQ) > diam(T^) — 1. 

Proof. The existence of an ancestor Tq comes from Lemma f3.2[ so we just have to show 
that diam(rp) > diam(r*) — 1. Tq is surrounded by a connected set of hexagons that 
touch Tq and whose spins are all the same (this is the "external boundary" of the set 
of hexagons that are in the interior of Tq). From this set, one can extract an m-path of 
constant sign whose diameter is bounded above by diaiu^) + 1. Since such a constant 
sign m-path is stable for the dynamics, T* must lie within its interior. This concludes the 
proof. □ 

At this point, we need to show that the above argument can be repeated and the 
construction done simultaneously for all curves in !Fg and J-g (for each n). First of all 
notice that, for a fixed e, any b-path T* of diameter less than e/2 can be approximated 
by a closest stable edge, provided that one is found within the ball of radius e/2 that 
contains the F*, with the probability of this last event clearly going to 1 as 5 — > 0, when 
we restrict attention to B l (6/e). For a b-path outside B 1 (6/e), we already noticed that 
it can be approximated by any other b-path also outside B l (6/e). As for the remaining 
b-paths, notice that the total number of boundaries that intersect the ball B l (6/e) cannot 
exceed const x (s5)~ 2 (in fact, the total number of pieces in which the boundaries that 
intersect B 1 (Q/e) are divided cannot exceed const x (e5)~ 2 ). So, we can carry out the 
above construction simultaneously for all the boundaries that touch B 1 (6/e), having to 
deal with at most const x (e5)~ 2 segments of b-paths of diameter of order at least e. 
Therefore, letting Y 5 n = dist^^, J~s), we can apply once again Lemma and conclude 
that 

H°(Y? >e)< {e5)- 2 e- c " {£/&) . (19) 

To show that Y s n — > as 5 —>■ /x°-almost surely and thus conclude the proof, it 
suffices to show that, Ve > 0, //(limsup^Q Y s n > e) = 0. To that end, first take a 
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sequence 5k = l/2 fe and notice that 

£/A^>-)<£^-^<oc, ( 2 °) 

k=0 k=0 

where we have made use of (fTTIJ) . Equation (ffllj) implies that we can apply the Borel- 
Cantelli lemma and deduce that ^(Hmsup^^ 17* > e) = 0, We > 0. In order to handle 
the values of 5 not in the sequence 5k, that is for those 5 such that 5k+i < 5 < 5k for some 
k, we use the following double bound, valid for any < a < 1, 

ad(u, v) < d(au, av) < — d(u, v), (21) 

a 

which implies that aY^ < Y£s k — q^"' ^ ne ^ wo bounds in equation (J21|) come from 
writing d(cra, av) as d(au, av) = inf^/ |0'| 2 ) -1 ds 1 = a inf^, + a 2 |0| 2 ) -1 ds, where 

4>'{s') are smooth curves joining au with av, while <fi(s) are smooth curves joining u with 
v. 

The proof of the theorem is now complete. □ 



3.2 Proof of Corollary I2TT1 

The Corollary is an immediate consequence of Theorem ^ and of the following general 
fact, of which we include the proof for completeness. 

Lemma 3.5. If {X$\, {Yg} (f or 5 > 0), and X are random variables taking values in a 
complete, separable metric space S (whose a -algebra is the Borel algebra) with {X$} and 
{Ys} all defined on the same probability space, then if Xg converges in distribution to X 
and the metric distance between X$ and Ys tends to zero almost surely as 5 — > ; Y$ also 
converges in distribution to X . 

Proof. Since X$ converges to X in distribution, the family {X$} is relatively compact 
and therefore tight by an application of Prohorov's Theorem (using the fact that S is 
a complete, separable metric space - see, e.g., [7j). Then, for any bounded, continu- 
ous, real function / on S, and for any e > 0, there exists a compact set K such that 
f \f (Xs)\I{x s <£K}dP < £ and f \f (Yg)\I{x s ^K}dP < £ for all 5, where /{.} is the indicator 
function and P the probability measure of the probability space of {X$} and {Ys}. Thus, 
for small enough 5, 

| J f(X 5 )dP- J f(Y s )dP\ < J \f(X s ) - f(Y 5 )\I {XseK} dP + 2e < 3e, (22) 

where in the last inequality we use the uniform continuity of / when restricted to the 
compact set K and the fact that the metric distance between Xs and Ys goes to as 
*-0. n 

To conclude the proof of the Corollary, it is enough to apply Lemma 13.51 to {^}s, 
{lJ^}s, Hsi ( or 5 to be more precise, to the random variables of which those are the dis- 
tributions), for each n G [l,oo], where fi s i is the unique scaling limit of critical site 
percolation on the triangular lattice, p 
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4 Dependent site percolation models on the hexago- 
nal and triangular lattice 

The model that we have presented and discussed in Section El has been chosen as a sort 
of paradigm, but is not the only one for which such results can be proved. In fact, it is 
not the original model for which such results were obtained. 

In this section we describe some percolation models on the hexagonal lattice and 
prove that they have the same scaling limit as critical (independent) site percolation on 
the triangular lattice. None are independent percolation models, but nonetheless, they 
represent explicit examples of critical percolation models on different lattices with the 
same scaling limit. Besides, the construction of the models on the hexagonal lattice can 
be seen as a simple and natural way of producing percolation models for which all the 
sites of the external (site) boundary of any constant sign cluster C belong to a unique 
cluster C of opposite sign. In other words, this implies that the the boundaries between 
clusters of opposite sign form a nested collection of loops, a property that site percolation 
on the triangular lattice possesses automatically because of the self-matching property of 
T (which is crucial in Smirnov's proof of the existence and uniqueness of the scaling limit 
of crossing probabilities and Cardy's formula). 

4.1 The models 

The percolation models that we briefly describe here can be constructed by means of 
a natural zero-temperature Glauber dynamics which is the zero-temperature case of Do- 
many 's stochastic Ising ferromagnet on the hexagonal lattice ^3] • The cellular automaton 
(i.e., Domany's stochastic Ising ferromagnet at zero temperature) that gives rise to those 
percolation models can also be realized on the triangular lattice with flips when a site 
disagrees with six, five and sometimes four of its six neighbors. The initial state a con- 
sists of an assignment of —1 or +1 with equal probability to each site of the hexagonal or 
triangular lattice (depending on which version of the cellular automaton we are referring 
to). In the first version, H, as a bipartite graph, is partitioned into two subsets A and 
B which are alternately updated so that each o~ x is forced to agree with a majority of 
its three neighbors (which are in the other subset). In the second version, all sites are 
updated simultaneously according to a rule based on a deterministic pairing of the six 
neighbors of every site into three pairs. The rule is that a x flips if and only if it disagrees 
with both sites in two or more of its three neighbor pairs; thus there is (resp., is not) 
a flip if the number D x of disagreeing neighbors is > 5 (resp., < 3) and there is also a 
flip for some cases of D x = 4. These percolation models on EI and T are investigated 
in ^0], where Cardy's formula for rectangular crossing probabilities is proved to hold in 
the scaling limit. 

The discrete time cellular automaton corresponding to the zero-temperature case of 
Domany's stochastic Ising ferromagnet on the hexagonal lattice can be considered as a 
simplified version of a continuous time Markov process where an independent (rate 1) Pois- 
son clock is assigned to each site i6l, and the spin at site x is updated (with the same 
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rule as in our discrete time dynamics) when the corresponding clock rings. (In particular, 
they have the same stable configurations - see Figure EJ) The percolation properties of the 
final state <t°° of that process were studied, both rigorously and numerically, in [22]; the 
results there (about critical exponents rather than the continuum scaling limit) strongly 
suggest that that dependent percolation model is also in the same universality class as 
independent percolation. Similar stochastic processes on different types of lattices have 
been studied in various papers. See, for example, [HI EH U3 E21 E31 EH EH] for models 
on Z d and [21] for a model on the homogeneous tree of degree three. Such models are 
also discussed extensively in the physics literature, usually on 7L d (see, for example, [TH] 
and |30j). Numerical simulations have been done by Nienhius [HE] and rigorous results for 
both the continuous and discrete dynamics have been obtained in [U|, including a detailed 
analysis of the discrete time (synchronous) case. 

Let us now describe in more detail the two deterministic cellular automata (on EI and 
T). Later, we will prove the equivalence of the percolation models generated by those 
cellular automata. 

Zero-temperature Domany model 

Consider the homogeneous ferromagnet on the hexagonal lattice HI with states denoted 
by o = {cr x }xew, &x = ±1, and with (formal) Hamiltonian 

where ^2{x y) denotes the sum over all pairs of neighbor sites, each pair counted once. We 
write M m {x) for the set of three neighbors of x, and indicate with 

A x H(a) = 2 a * a v ( 24 ) 

the change in the Hamiltonian when the spin a x at site x is flipped (i.e., changes sign). 

The hexagonal lattice HI is partitioned into two subsets A and B in such a way that 
all three neighbors of a site x in A (resp., B) are in B (resp., A). By joining two sites 
of A whenever they are next-nearest neighbors in the hexagonal lattice (two steps away 
from each other), we get a triangular lattice (the same with B) (see figure |7J). The 
synchronous dynamics is such that all the sites in the sublattice A (resp., B) are updated 
simultaneously. 

We now define the discrete time Markov process a n , n G N, with state space Sm = 
{ — 1, +1} H , which is the zero temperature limit of a model of Domany J^l, as follows: 

• The initial state a is chosen from a symmetric Bernoulli product measure. 

• At odd times n — 1,3,..., the spins in the sublattice A are updated according to 
the following rule: a x , x G A, is flipped if and only if A x H(a) < 0. 

• At even times n — 2, A, ... , the spins in the sublattice B are updated according to 
the same rule as for those of the sublattice A. 
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Figure 6: Example of a (local) stable configuration for the zero-temperature Domany 
dynamics. Heavy lines on edges of EI connect, say, plus spins, while heavy broken lines 
connect minus spins. The dotted lines drawn on edges of the dual lattice are the perpen- 
dicular bisectors of unsatisfied edges, indicating the boundaries between plus and minus 
clusters. Every spin has at least two neighbors of the same sign, since only loops and 
barbells are stable under the effect of the dynamics. 



Cellular automaton on T 

We define here a deterministic cellular automaton Q on the triangular lattice T, with 
random initial state chosen by assigning value +1 or —1 independently, with equal prob- 
ability, to each site of T. 

Given some site x G T, group its six T-neighbors y in three disjoint pairs 
{2/35 2/4 {vtiV%}i so that yf and y% are T-neighbors, and so on for the other two pairs. 
Translate this construction to all sites x G T, thus producing three pairs of sites {yf, y%}, 
{2/3 ) 2/f }> {2/5 > Vt } associated to each site x G T. (Note that this construction does not need 
to specify how T is embedded in M 2 .) Site x is updated at times m — 1, 2, . . . according 
to the following rule: the spin at site x is changed from a x to — a x if and only if at least 
two of its pairs of neighbors have the same sign and this sign is — a x . 

4.2 Equivalence between the models on EI and T 

We show here how the models on the hexagonal and on the triangular lattice are related 
through a star-triangle transformation. More precisely, we will show that the dynamics 
on the triangular lattice T is equivalent to the alternating sublattice dynamics on the 
hexagonal lattice EI when restricted to the sublattice B for even times n = 2m. 

To see this, start with T and construct an hexagonal lattice EI' by means of a star- 
triangle transformation (see, for example, p. 335 of [TH]) such that a site is added at 
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the center of each of the triangles (x, yf, yf), (x, y%, y%), and {x^yf,y^) (the sites yf are 
defined in the previous subsection). HI' may be partitioned into two triangular sublattices 
A 1 and B' with B' = T. One can now see that the dynamics on T for m — 1,2,... and 
the alternating sublattice dynamics on EI' restricted to B' for even times n = 2m are the 
same. 




Figure 7: A star-triangle transformation. 



An immediate consequence of this equivalence between the two cellular automata is 
that the two families of percolation models that they produce are also equivalent in an 
obvious way through a star-triangle transformation. To be more precise, the percolation 
models defined on T by Q for times m = 1, 2, . . . are the same as those defined on B by 
the zero-temperature Domany model for even times n = 2m. 

4.3 Results for the zero-temperature Domany model 

In this section we explain how the results obtained for the percolation models /i^ gen- 
erated by the cellular automaton T are also valid for the percolation models defined by 
the zero-temperature Domany model. Because of the results of the previous section, we 
can actually consider the percolation models jx^ on ST generated by Q, for which we 
have results analogous to Theorem ^ and Corollary 12.11 Results equivalent to Propo- 
sitions 12.11 12.21 and 12.31 are contained in The proof of the main theorem (i.e., the 
analogue of Theorem ^) is basically the same as for the models generated by T, so we just 
point out the differences. We follow here the setup and notation of P, but give all the 
relevant definitions in order to make this paper self-consistent. 

Let us consider a loop T in the triangular sublattice B, written as an ordered sequence 
of sites (y , yi, . . . , y^) with k > 3, which are distinct except that y^ = y . For i — 1, . . . , k, 
let Q be the unique site in A that is an H-neighbor of both y^i and y^. We call T an 
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s-loop if (i, . . . , (k are all distinct. Similarly, a (site-self avoiding) path (yo, yi, . . . , yu) in 
B, between y and y^, is called an s-path if Ci> • • & are an distinct. Notice that any 
path in B between y and y' (seen as a collection of sites) contains an s-path between y 
and y' . An s-loop of constant sign is stable for the dynamics since at the next update of 
A the presence of the constant sign s-loop in B will produce a stable loop of that sign 
in the hexagonal lattice. Similarly an s-path of constant sign between y and y' will be 
stable if y and y' are stable — e.g., if they each belong to an s-loop. A triangular loop 
Xi,X2,x% G B with a common H-neighbor £ G A is called a star; it is not an s-loop. A 
triangular loop in B that is not a star is an s-loop and will be called an antistar, while 
any loop in B that contains more than three sites contains an s-loop. 

With these definitions, the proof of the main theorem for the zero-temperature Domany 
model is the same as that presented in Section 01 for our prototypical model, with the role 
of loops in the original proof played here by s-loops (in particular antistars), and that of 
m-paths by s-paths (see Section |3] below). 

4.4 An amusing further example: totally synchronous dynamics 
on H 

The model that we consider here corresponds to the zero-temperature Glauber dynamics 
on the hexagonal lattice with all sites updated simultaneously at discrete times, without 
alternating between two subsets (contrary to the case of the zero-temperature Domany 
model), with initial configuration a chosen according to a symmetric Bernoulli product 
measure. 

Let us partition EI in two subsets A and B as before and define the family of percolation 
models {ft™, m = 0, 1, . . .}, where p,™ is the distribution of a 2m (at even times) restricted 
to the subset A (naturally endowed with a triangular lattice structure, so that Ji° A is the 
distribution of critical site percolation). 

The main difference consists in the fact that a n does not fixate as n — > oo since there is 
a positive density of spins that flip infinitely many times. To see this, consider a loop V in 
EI containing an even number of sites and such that at time zero its spins are alternately 
plus and minus. At any time n, every spin in T has two neighbors of opposite sign and 
will therefore flip at the next update. Thus, the spins in T never stop flipping. 

However, there is a simple observation that tremendously simplifies the analysis of 
this totally synchronous dynamics. Namely, that if we restrict attention to sublattice A 
at odd times n = 1, 3, 5, . . . and sublattice B at even times n — 0, 2, 4, . . ., the dynamics 
is identical to the zero-temperature Domany dynamics discussed above (let us call this 
a"). On the other hand, if we instead observe A at even times and B at odd times, this is 
identical to an alternative zero-temperature Domany type dynamics cx^ but with the first 
(and third and . . .) update on B rather than A. Furthermore, (cr2)^=o an d (ffeOn^o are 
completely independent of each other. We conclude that there are two distinct limits cr^° 
and (independent of each other) and that the scaling limit of a n restricted to either 
A or B is the same as for independent critical percolation on T. But it appears that for 
any n, a n on all of EI should be subcritical and thus have a trivial scaling limit. 
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